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It is commonly believed that the fidelity of quantum teleportation using localized quantum objects 
with one party or both accelerated in vacuum would be degraded due to the heat up by the Unruh 
effect. In this paper we point out that the Unruh effect is not the whole story in accounting for 
all the relativistic effects in quantum teleportation. First, there could be degradation of fidelity 
by a common field environment even when both quantum objects are in inertial motion. Second, 
relativistic effects entering the description of the dynamics such as frame dependence, time dilation, 
and Doppler shift, already existent in inertial motion, can compete with or even overwhelm the 
effect due to uniform acceleration in a quantum field. We show it is not true that larger acceleration 
of an object would necessarily lead to a faster degradation of fidelity. These claims are based on 
four cases of quantum teleportation we studied using two Unruh-DeWitt detectors coupled via a 
common quantum field initially in the Minkowski vacuum. We find the quantum entanglement 
evaluated around the light cone, rather than the conventional ones evaluated on the Minkowski 
time slices, is the necessary condition for the averaged fidelity of quantum teleportation beating 
the classical one. These results are useful as a guide to making judicious choices of states and 
parameter ranges and estimation of the efficiency of quantum teleportation in relativistic quantum 
systems under environmental influences. 

PACS numbers: 04.62,+v, 03.67.-a, 03.65.Ud, 03.65.Yz 


I. INTRODUCTION 

Quantum teleportation (QT) is by now quite well recognized as a feature process in the application of quantum 
information w- a novel and exclusively quantum process QT is also of basic theoretical interest because it necessi¬ 
tates a proper treatment of quantum measurement and entanglement dynamics in realistic physical conditions, such 
as environmental influences. The advent of a new era of quantum sciences and engineering demands more precise 
understanding and further clarification of such fundamental issues. This includes quantum information and classical 
information, quantum nonlocality and relativistic locality, and spacelike correlations and causality. The study of these 
issues in a relativistic setting now belongs to a new field called relativistic quantum information [3]. 

The first scheme of QT was proposed by Bennett et al. (BBCJPW) [5], in which an unknown state of a qubit 
C is teleported from one spatially localized agent Alice to another agent Bob using an entangled pair of qubits A 
and B prepared in one of the Bell states and shared by Alice and Bob, respectively. Such an idea was then adapted 
to the systems with continuous variables such as harmonic oscillators (HOs) by Vaidman [S], who introduced an 
ideal Einstein-Podolsky-Rosen (EPR) state for the shared entangled pair to teleport an unknown coherent state. 
Braunstein and Kimble (BK) [5] generalized Vaidman’s scheme from the ideal EPR states with exact correlations to 
squeezed coherent states. In doing so the uncertainty of the measurable quantities has to be considered, which reduces 
the degree of entanglement of the AB pair as well as the fidelity of quantum teleportation (FiQT). 

Alsing and Milburn made the first attempt of calculating the FiQT between two moving cavities in relativistic 
motions [9]-one is at rest (Alice), the other is uniformly accelerated (Bob, called Rob in [5] with the initial “R” for 
“Rindler observer”; we follow this convention in Section [V] for a similar setup) in the Minkowski frame-to see how the 
fidelity is degraded by the Unruh effect (also see Refs. [101 ITT] ). Later Landulfo and Matsas considered a complete 
BBCJPW QT in a two-level detector qubit model, where Rob’s detector is uniformly accelerated and interacting with 
the quantum field only in a finite duration. They found that the FiQT in the future asymptotic region using the out 


* Electronic address: sylin@cc.ncue.edu.tw 
tElectronic address: chouch@mail.ncku.edu.tw 
^ Electronic address: blhu@umd.edu 



2 


state of the entangled pair is indeed reduced by the Unruh effect experienced by Rob [T2]. Along the line of Ref. [5] 
Friis et al. [T3] studied the role of the dynamical Casimir effect in the QT between cavities in relativistic motions. 

Alternatively, Shiokawa P3] considered QT in the Unruh- DeWitt (UD) detector theory [151 [TB] with the agents in 
motions similar to those in Ref. [5] , but based on the BK scheme in the interaction region : An unknown coherent state 
of a UD detector with internal HO is teleported from Alice to Bob using an entangled pair of similar UD detectors 
initially in a two-mode squeezed state and shared by Alice and Bob. Unfortunately, the FiQT considered in Ref. [14] 
is not the physical one. More careful consideration is needed to get the correct results PH- 

Indeed, when considering QT in a fully relativistic system, particularly in the interaction region of the localized 
objects and quantum fields, one has to take all the factors listed below into account consistently. 


A. Relativistic effects 

Localized objects in a relativistic system may behave differently when observed in different reference frames: 

a. Frame dependence Since quantum entanglement between two spatially localized degrees of freedom is a kind of 
spacelike correlation in a quantum state, which depends on reference frames, quantum entanglement of two localized 
objects separated in space is frame dependent. 

b. Time dilation When two localized objects in uniform motion have a nonzero relative speed, both will perceive 
the same time dilation of each other in their rest frame constructed by the radar times and distances. If one object 
undergoes some phase of acceleration but the other does not, e.g., the worldlines in the twin problem, then the time 
dilations perceived by these two objects will be asymmetric. All these time dilation effects are included in the proper 
time parametrization of the worldline of an object localized in space. 

c. Relativistic Doppler shift Suppose Alice continuously sends a clock signal periodic in her proper time to Bob, 
then Bob will see Alice’s clock running slower or faster than the one at rest when the received signal is redshifted or 
blueshifted, depending on their relative motion. 

These three basic properties of relativistic quantum systems essential for the consideration of QT have not been 
properly recognized or explored in detail or depth. 


B. Environmental influences 

The qubits or detectors in question are unavoidably coupled with quantum fields, which act as an ubiquitous 
environment: 

d. Quantum decoherence Each qubit or HO can be decohered by virtue of its coupling to a quantum field. 
However, mutual influences mediated by the field between two localized qubits or HOs when placed in close range 
can lessen the decoherence on each. 

e. Entanglement dynamics The entanglement between two qubits or HOs changes in time as their reduced state 
evolves. 

/. Unruh effect A pointlike object such as a UD detector coupled with a quantum field and uniformly accelerated 
in the Minkowski vacuum of the field would experience a thermal bath of the field quanta at the Unruh temperature 
proportional to its proper acceleration m- 


C. New issues in dealing with quantum teleportation 

The above factors have been considered earlier in some detail in our study on entanglement dynamics [18H20j . but 
there are new issues of foundational value that need be included in the consideration of QT. Below we mention three 
issues related to relativistic open quantum systems: 

g. Measurement in different frames Quantum states make sense only in a given frame in which a Hamiltonian 
is well defined [SUES]- Two quantum states of the same system with quantum fields in different frames are directly 
comparable only on those totally overlapping time slices associated with some moment in each frame. By a measure¬ 
ment local in space, e.g., on a pointlike UD detector coupled with a quantum field, quantum states of the combined 
system in different frames can be interpreted as if they collapsed on different time slices passing through the same 
measurement event [41] . Nevertheless, the postmeasurement states will evolve to the same state up to a coordinate 
transformation when they are compared at some time slice in the future, if the combined system respects relativistic 


covariance 
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h. Consistency of entangled pair As indicated in the BK scheme, the FiQT could depend on (i) quantum entan¬ 
glement of the entangled pair and (ii) the consistency of the quantum state of the entangled pair with their initial 
state. Both would be reduced by the coupling with an environment, and applying an improved protocol of QT may 
suppress the deflection of (ii). 

i. Comparing FiQT and entanglement It is easy to modify the BBCJPW scheme to see that the FiQT of qubits 
in pure states can be either 1 or 0, depending only on whether the qubit pair is entangled or not. In contrast to 
qubits in pure states, the best possible FiQT in the BK scheme depends on how strong the HO pair is entangled [24| . 
To compare the degree of entanglement of the entangled pair and the FiQT applying them in relativistic systems, 
Shiokawa considered a “pseudofidelity” of QT evaluated on the same time slice for the degree of entanglement by 
imagining that right at the moment Alice has just performed the joint measurement, Bob gets the information of the 
outcome from Alice instantaneously and immediately performs the proper local operations on his part of the entangled 
pair mm- In reality, classical information needs some time to travel from Alice to Bob, and during the traveling 
time, Bob’s part of the entangled pair keeps evolving, so the physical FiQT will not be equal to the “pseudofidelity” and 
is thus incommensurate in general with the degree of entanglement of the entangled pair evaluated on the Minkowski 
time slice. This feature has been overlooked in the literature. 


D. Organization of this paper 

To address all the above issues consistently and thoroughly, we start with the action of a fully relativistic system. 
We introduce the model in Sec. II, then derive the formula of the FiQT for our model in Sec. Ill, where we discuss 
the relation between the fidelity and the degree of quantum entanglement of the detector pair. In Secs. IV to VII, 
respectively we apply our formulation to four representative cases with Alice at rest and 1) Bob also at rest |18| , 
2) Bob (Rob) uniformly accelerated in a finite period of time O |T5] > 3) Bob being the traveling twin in the twin 
problem [25], and 4) Bob undergoing alternating uniform acceleration |26j . The trajectories and kinematics of each 
case can be found in the sample references given above. Finally we summarize and discuss our findings in Sec. VIII. 
In the Appendix we show the consistency of the reduced states of the detectors under the spatially local projective 
measurements. 


II. MODEL 


Consider a model with three identical Unruh-DeWitt detectors A , B , and C moving in a quantum field 4>(:r) in 
(3+l)-dimensional Minkowski space. The internal degrees of freedom Q A , Q B , and Q c of the pointlike detectors A, 
B , and C, respectively, behave like simple harmonic oscillators with mass m = 1 and natural frequency Cl. The action 
of the combined system is given by m 


1 r 


S = - I d i x^g-d l _ l <^(x)d^$(x) + ^ / dr d - ( d d Q d ) 2 - C%Q\ 


d =A,B,C" 


+ ^2 A 0 f d 4 x f dr d Q d (T d )$(x)5 4 (x M - z^(r d )) 
d=A,B 


(1) 


where p = 0, 1,2,3; g^p = diag(— 1, 1, 1,1); d d = d/dr A ; r A , r B and r G are proper times for Qa , Qb , and Qc, 
respectively; and the lightspeed c = 1. The scalar field 4> x (f) = 4>(t,x) = 4>(;r) is assumed to be massless, and Ao is 
the coupling constant. Detectors A and B are held by Alice and Bob, respectively, who may be moving in different 
ways, while detector C carries the quantum state to be teleported and goes with the sender. 

Suppose the initial state of the combined system defined on the t = 0 hypersurface in the Minkowski coordinates is 
a product state ® p AB ®> p^’ r °\ where p $ = 0m) (0m| is the Minkowski vacuum of the field, p AB is a two-mode 

squeezed state of detectors A and B, and p^’^ is a squeezed coherent state of detector C with a = a R + ia : and 
r 0 the squeezed parameter. In the (A', A) representation [231 Ell (the double Fourier transform of the usual Wigner 
function, namely the “Wigner characteristic function” [28]), we express the last two as 


p ( b a ’ r °\K c ,A c ) = J dH c e^ C (Q c \p^ ro) \Q' c ) 


= exp 


A ( A e *»(i < c f + W 2 '"(a c ) 2 ) + i ( \I^o r k° - Vmta.A? 


Cl 


2 h V 2 Cl 


q c ,q' c =y,ct( AC/2) 
2h 


( 2 ) 
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and 


Pab (K a ,K b ,A a ,A b ) 


1 

= exp — — 


8 IP 


1 


or 


~^(K A + K B f + ^(A A + A B f + ^ ( K A - K B f + A_(A a - A b ) 2 


a* 


(3) 


with parameters a and /3. One may choose a = e ri y/U/Sl and (3 = e ri y/hCl, where rq is the squeezed parameter. 
As n — > oo, p AB goes to an ideal EPR state with the correlations (Qa — Qb) = (Pa + Pb) =0 without uncertainty, 
while Q a + Qb and P A — P B are totally uncertain. Here, Pd is the conjugate momentum to Q d . 

In general the factors in pjP’ r °\K c ,A C ) will vary in time. To concentrate on the best FiQT that the entangled 
AB pair can offer, however, we follow Ref. |14| and assume the dynamics of is frozen or, equivalently, assume 

,r °^ is created just before teleportation. 

At t = 0 in the Minkowski frame, the detectors A and B start to couple with the field, while the detector C 
is isolated from others. By virtue of the linearity of the combined system Q, the quantum state of the combined 
system started with a Gaussian state will always be Gaussian, and therefore the reduced state of the three detectors 
is Gaussian for all times. In the (AT, A) representation the reduced Wigner function at the coordinate time x° = T in 
the reference frame of some observer has the form 


Pabc( k < a ; t ) = exp 


d (T))K d - (P d (T))A d 


d,d' 


where d, d' = A, B,C , and the factors 


(k a Q aa ,(T)K a ' + A d Pdd'(T)A d ' - 2A' d Pdd'(T)A d ' 


Qdd'(T) 

Pdd'(T) 

ftdd'Cn 


ns ns 


i6K d i8K d ' 
ihS ihS 


Pabc 


SA d SA d 

ns ms 

iSK A dA d 


r Pabc 

tPabg 


= (SQ d (T d (T)),SQd'(Td'(T))}, 

K=A—0 

= <6Pd(Td(T)),5Pd>(Td'(T))), 


K=A=0 


= (6Qd(T d (T)),SP d ,(T d ,(T))), 


K=A=0 


(4) 

(5) 

( 6 ) 
(7) 


are actually those symmetric two-point correlators of the detectors in their covariance matrices ((O^O') = (OO' + 
O'0)/2 and SO = O — (O)), which can be obtained in the Heisenberg picture by taking the expectation values of the 
evolving operators with respect to the initial state defined on the fiducial time slice. 


III. FIDELITY OF QUANTUM TELEPORTATION AND ENTANGLEMENT 


For our later use, below we reexpress and generalize the definitions and calculations for QT of a Gaussian state 
from Alice to Bob in Refs. EH ED] in terms of the ( K , A) representation. Suppose the reduced state of the 
three detectors continuously evolves to p ABC (K, A ; p) in the Minkowski frame when Alice’s and Bob’s proper times 
are t a = r A (t\) and t b = r B (t\), respectively. At this moment Alice preforms a joint Gaussian measurement 
locally in space on A and C so that the postmeasurement state right after t\ in the Minkowski frame becomes 
p ABC (K, A; fi) = p A c(K A ,K°, A A ,A c )p B (K B , A B ), where we assume the quantum state of detectors A and C 
becomes another two-mode squeezed state 


Pac( kA > a ' C , A 4 , A c ) = exp 


mv^ c -^ aC 


1 


(k™ Q 


-—2 ( K m Q mn K n + A m V mn A n - 2K m U mn A' 


( 8 ) 


with m,n = A,C so that Alice gets the outcome (3 = (3r + i/3i. (Here and below, the Einstein notation of summing 
over repeated dummy indices is understood, and Y) m n ignored.) Then Eq. (l8|) yields the reduced state of detector 
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B 


~ (rB] _ at [d 2 K G <PlC A o,). A c 
p B (JC ) - N b J — —p^yc ,/C )x 


27 rh 27 tTi 
P ABC {K a ,K b ,K c -M), 


(9) 


right after t b , where N B is the normalization constant, /C d = (A' d ,A d ), and d 2 KL d = dI\ d dA d . If we require 
1 = Tr s p B (= P b \k b =a b =o)> then N B will depend on /3. Alternatively, following Ref. [Hj, we can require N B to be 
independent of /3, and then Tr B p B will be proportional to the probability P(/3) of finding detectors A and C in the 
state fel). Let Tr s p B = P(/3); then, the normalization condition reads [0 = (0,0)] 


1 = / d 2 ldP(P) = / dP R dfr p B (IC B = 0) 


= N b / d(3 B d/3i 
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= N b 


n b 
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exp 
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= / d 2 IC A exp —^ 
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2 h 


2 h z 


(q 1 1 ] a + Qaa) (K a ) 2 + (v [ l\ + V AA ) {A A Y - 2K a ( 


- 2K A (n l l ] A + h AA ) A 


after inserting Eqs. (|4]) and ^ into the integrand. Here, S I n l denotes the value of the factor S = Q, V, or 1Z being 
taken at t n — e with e —» 0+. Thus, we have 


Nb = ^ ^(q 1 a ] a + Qaa) (v^a + Paa) - (tL^ + Paa) 2 ■ ( 10 ) 

Right after the joint measurement on A and C, Alice sends the outcome /3 of the measurement to Bob by a 
classical signal at the speed of light. Suppose the signal reaches Bob at his proper time t b = T° dv = T adv (t\) 
[here, “adv” stands for “advanced” [31], and r adv is the advanced time defined by | z B {T adv (t)) — z A (t) | 2 = 0 with 
z° B (r adv {t )) > z A (t)\, when the reduced state of detector B has evolved from the postmeasurement state (|9]) to p' B . 
According to the information received, Bob could choose a suitable operation on detector B to turn its quantum state 
to a copy of the original unknown state carried by detector C. In the BK scheme mm, the operation Bob should 
perform is a displacement by /? in the phase space of detector B , namely, p out = D(/3)p" Bl where p' B is the reduced 
state of detector B keeps evolving from rf 6 * 1 ' to the operation event, and D(f 3) is the displacement operator, or in the 
{K, A) representation, 


Pout(IC B ) = p B (K. B ) ex p 


P r K b - v / 2Ml/3/A i 


( 11 ) 


The fidelity of quantum teleportation (FiQT) from |a,ro) c to |a,ro) B is then defined as 

= b( a , r 0 I Pout I r o) B 
B Pout 


( 12 ) 


If we have an ensemble of the distinguishable ABC triplets of the detectors, the quantity we are interested in will be 
the averaged FiQT [32], defined by 


F av = / d 2 (3P(f3)F((3) = / dfadfa 


Tr bPb 

Tr B ~p" B 


B ( a , r 0 \p ou t\a, ro) B , 


(13) 


since Tr B p out = p out {IC B = 0) = p B (IC B = 0) = Tr B p” B . 
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FIG. 1: Setup for QT from Alice (thick dashed worldline) to Bob (thick solid worldline); both are at rest in the Minskowski 
vacuum. The gray solid curve represents the t\ slice in some coordinate system, and the gray dashed horizontal lines represent 
the t slices in the Minkowski coordinates. The shaded region represents the future light cone of the joint measurement event 
on A and C by Alice (red cross). 


A. Direct comparison of FiQT and entanglement 


In Ref. 


Mari and Vitali showed that the optimal averaged FiQT of a coherent state is bounded above by 


F op t < 


1 


1 + (2c-/h) ’ 


(14) 


where c_ is the lowest symplectic eigenvalue of the partially transposed covariance matrix in the reduced state of the 
entangled AB pair defined on the time slice right before the joint measurement at t\ USES. c_ can be related to quan¬ 
tum entanglement of the AB pair by noting that the logarithmic negativity is given by Ej i/ = max{0, — log 2 (2c_/fi.)}. 
Nevertheless, the dynamics of detector B between Alice’s measurement and Bob’s operation have been ignored in 
obtaining the above inequality. In a relativistic open quantum system, (14) does not make exact sense, since the 


averaged FiQT on the left side of (14) is a timelike correlation connecting the joint measurement event by Alice and 


the operation event by Bob, while the quantity on the right side of (14) is a spacelike correlation extracted from the 
covariance matrix of detectors A and B defined on the hypersurface of simultaneity right before the wave functional 
collapses. 

To compare the averaged FiQT directly with a function of c_ defined on the ti-slice in the Minkowski frame, one 
might imagine that Bob receives the outcome /3 and makes the proper operation on detector B instantaneously at t b 
when the worldline of B intersects the fi-slice (see Fig. 0 El, which is unphysical. 

A better way to make a direct comparison is to transform the combined system to a new reference frame with the 
fiducial time slice overlapping with the t = 0 hypersurface in our original setup but the time slice passing Alice’s 
measurement event being very close to the future light cone of the event (e.g., the gray solid curve in Fig. [l] joining 
Alice’s worldline at rj 4 and Bob’s worldline at t[ b = T^ dv — e, e —>■ 0+). Then the wave functional defined in this 
new reference frame is collapsed around the future light cone of the joint measurement event, right after which Bob 
receives the signal from Alice and immediately performs the operation on detector B (at T? dv + e in Fig. 0, which is 


still around the same future light cone and so p" B ss p' B ~ p B - I n this way both sides of (|14j) are evaluated around the 
future light cone of Alice’s measurement event, or around the past light cone of Bob’s operation event, and both sides 
of (14) will be independent of the reference frame in a relativistic detector-held system when e —» 0. In the Appendix 


we show that the reduced state of detector B collapsed around the light cone of the joint measurement event on A and 
C is consistent with the reduced state initiated with the one collapsed simultaneously with the measurement event in 
a conventional reference frame and then evolves to the future light cone of the event. Actually, the reduced state of 
detector B at the moment that Bob is crossing the future light cone of Alice’s spatially local measurement event is 
independent of the choice of coordinates here. 

Denoting by t' the coordinate time of a new coordinate system such that r j4 (t , 1 ) = r A (ti) = rj 4 and T B (t[) = 


,—adv 

T 1 


— e = 


= t[ b at t' = t[, and assuming f/ B {r^ dv + e) ~ P B { T \ dv - e) as e -> 0+. Then, we can repeat the same 
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approach described earlier in this section to reduce Eq. (131 to 


F av = J d/3 R d/3i B (a,r 0 \pout\a,r 0 } B = J d/3 R d/3 I C ^^-p ( g ,ro) * (K. B )p out (IC B ) 1 


(15) 


where p^ Ao ^ in the (K, A) representation is the same as Eq. ([ 2 1 except the index C there is replaced by B. From 
Eqs. Q and (11), with the help of Eqs. @ and ([8]), and with t\ replaced by t\ , we have 

Fav = N b J dpRdpj P abc (K, A; t[) x 


- * B ) - - Pi)(A c - A b ) 

A_ e 2r 0 (K B f + ^ne~ 2r °(A B ) 2 + K m Q mn K n + A m V mn A n - 2K m K mn A r ‘ 


= Nb J ^^Pabc(K,A; t[)(2n) 2 S 



(K c -K B )\ 5 



(A c - A B ) x 


6XP {'i 


A_ e 2r 0 (K B ) 2 + T ln-e 2r °(A B ) 2 + K m QmnK n + A m V mn A n - 2K m h mn A r ' 


(16) 


Thus, 


p — 
av 


h^nNs 

y/detV 


(17) 


where N B is the same as Eq. (10) except ti is replaced by t\ and 


( Q^ 1 ' 

/ ^ A A 


Qaa 
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F^aa FaA 
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F^aa — Fa a 


’ i n- 

>[i' 


pM 
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Vaa 


F B a Fqa 
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*BA 

[I'l 
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F^ba Fqa 
Q cc + he 2r °tt 
-F [ Pb-Fcc 


-1 


F^ab Fac 
•pM 

1 AB 

-npB-Fcc 


Vac 


(18) 


v [ B^ + Vcc + he- 2r °n ) 


Here the symmetric two-point correlators of the detectors, e.g., = Qdd'W) = {5Qd(T d (t , 1 ))6Qd'(T d '(t[)) are 

the expectation values of the operators of detector A at t a and the operators of detector B at t[ b = — e, with 

respect to the initial state of the combined system defined on the fiducial time slice t! = t = 0. One can easily write 
down a similar formula for the QT from Bob to Alice by switching their roles and letting detector C go with Bob. 

Note that F av in Eq. (17) is independent of a only if pp is in the form of Eq. (Jsj) , where the j3 terms are 
independent of K A or A A . The state Q is chosen so that the analytic calculation is the simplest while the result is 
still interesting. One may choose another state consistent with the ideal EPR state as the squeeze parameter r 2 —> 00 
instead; for example, K c and A c are replaced by ( K c — K A ) and (A c + A A ), respectively. Then, N B and the F av 
will be more complicated and will depend on a. In practice, the choice of the state may depend on the experimental 
setting. 

Below we consider the cases with the factors in the two-mode squeezed state ([8]) of detectors A and C right 
after the joint measurement given by Qaa = Qcc = ^ cosh2r 2 , Qac = ^sinh2r 2 , Vaa = Vcc = |Hcosh2r 2 , 
Vac = — |Osinh2r 2 with squeezed parameter r 2 , and lZ mn = 0. 

If the joint measurement on detectors A and C is done perfectly such that r 2 —> 00 , then from Eqs. ([17]), and 


(10), we have 


F av (rf,r[ B ) 


(he 2ro n~ l + (SQl >) ( he~ 2r °n + (SP 2 )\ - (< SQ-,6P+)y 


1-1/2 


(19) 


where Q- = Qa{t a ) — Qb{t[ b ) and P + = Pa{t a ) + Pb(t[ b ). If, in addition, the initial state p AB of detectors A 
and B in Eq. ([ 3 ]) were frozen in time and decoupled from the field, then one would have 


Fav(P,r ,B ) = F av ( 0,0) = 


sj (e 2r ° + e ~ 2ri ) (e~ 2r ° + e~ 2ri ) ’ 


( 20 ) 


























which implies that F av —► 1 as r\ —> oo when p AB is nearly an ideal EPR state, while F av —* 1/2 for ro = 0 as 
7 "i —>• 0 when p AB is almost the coherent state of free detectors. In the latter case F av = F c i = 1/2 is known as the 
best fidelity of “classical” teleportation of a coherent state carried by detector C using the coherent state of the AB 
pair 0 without considering the environmental influences. This does not imply that F av of QT must be greater than 
1/2, though. Once the correlations such as (Q-) = 0 needed in the protocol of QT become more uncertain than the 
minimum quantum uncertainty, F av — F c i will become negative. 

The degrees of quantum entanglement of the AB pair in their reduced state defined on the t\ slice, such as the 
logarithmic negativity E A r. can be evaluated by inserting the expressions for the two-point correlators of detectors A 
and B on that slice into the conventional formula nm eh ssj. Those correlators measure the correlations between 
the operators of detector A at some event (in Alice’s worldline at rf 4 ) and the operators of detector B at another 
event almost lightlike but still spacelike separated with the former (in Bob’s worldline at t[ b ). We call the quantum 
entanglement evaluated in this way as the “entanglement around the light cone” (EnLC). While the degrees of 
entanglement of two detectors obtained in the conventional ways depend on the choice of reference frames HE those 
for the EnLC do not. The inequality (14) implies that the EnLC between A and B (c_ < h/2 or i/y > 0) is a 
necessary condition for the averaged FiQT of coherent states to be better than the classical ones ( F opt > F c i). 


B. Ultraweak coupling limit 


In the ultraweak coupling limit, 7 = Aq/ 87 t is so small that 7 A 1 <C a, 17, where Ai corresponds to the time resolution 
or the frequency cutoff of our model p3j. From Eqs. (28)-(29), (32)-(33), and (B2)-(B8) in Ref. jT3| with a 2 = 
(h/V,)e~ 2ri and /3 2 = UVLe~ 2ri there (denoted by a and (3 in this paper), with 1 (qAi/fl) (7/17) O ( 7 A 1 /H) 2 , 
the elements of the covariance matrix for the AB pair at t\ with the initial state ©> can be approximated by 


gi 1 ') 

^■AA 

nd'l 


hC, _ 2 o"r. A 
217 6 


27T " + ((6Q A (Tf))X, v [ p « -nc ie - 2 ^ + ((SP A (X))X, 


^e- 2 ^ + (( 6Q b (t [ b )) 2 ) v , pM « hW ie -»ri B + (( 5P b {t' b ))X , 


Q l P « ^e-y^+^co S Q(X + T[ B ), 

[I'l 


217 


TV 


A B 


n BA ~ ^ S x e Arf+A B ) sin 17 (t / 4 + t' b ), 


n 


[i'l 

AA 


[I'l 

'BB 


( 21 ) 

( 22 ) 

(23) 

(24) 


up to h • 0(7/17). Here, C n = cosh2r n , S n = sinh2r„, ((<5Pj(r- 7 )) 2 ) v ss 17 2 ((^Q,-(r J ')) 2 ) v + v with j = A,B, and 
v = 2h'yAi/'rr. For simplicity, let us consider the cases with tq = 0 here. Then, Eq. (18) becomes 


V = 


hAX) 

0 


p{X,r[ B ) 

\ 

0 

|1 lA{rf) + v 


-|17*(t 1 V( b ) 

£>X(Tl A ,Ti B ) 



0 


§y(r? a b ) 

-fl7A(r( 4 ,r{ B ) 

0 

§17 B(t' b ) + v 

) 


h 4 0( 7/17), 


(25) 


where 


A(rf) = C 2 + e- 2 ^C 1 + 2im- 1 ((5Q A (X)) 2 ) Vl 
B{t’ b ) = 2 + 0 2 + e- 2 ^ B 0 1 + 217?r 1 (( ( 5Q B (r( B )) 2 ) v , 
X(X,t[ b ) = S 2 + e _7(riA+ ' riS ) cos!7(t ] a + t' b ) Si, 
y(X,T[ B ) = e- 7(T i A+r ^ ) sinl7(r 1 A + r( B )5i. 

So the averaged fidelity in the ultraweak coupling limit can be written in a simple form: 

2 A 


FaviX, r[ B ) = 


AB^{X 2 +y 2 ) 


0( 7 Ai/17). 


(26) 

(27) 

(28) 
(29) 


(30) 


Usually, (( SQj) 2 (r)) v ~ (±e 27T + constant) evolve smoothly in this limit, while 

X 2 + y 2 = Si + Sf e ~ 2 ^ + < B) + 2S\S 2 cos 17 (rj 4 + t' b ) 


( 31 ) 
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is oscillating in t 4 +t{ b due to the natural squeeze-antisqueeze oscillation of the two-mode squeezed state of detectors A 
and B at frequency f 1 m- The maximum (minimum) values of F av , denoted by F+ v (F av ), occur at cos « 1 

(— 1 ), when y = 0 and 


tj± f^A UB\ 
^av ( T 1 ^ ^1 ) 


2A 


AB- 


S 2 ± Si e~^ T i + < B '> 


(32) 


We call F+, the best averaged FiQT from Alice to Bob. 


C. Improved protocol 


Similar to the function of the local oscillators in the optical experiments of QT, if we perform a counter-rotation 
to p B in the phase space of ( Qb,Pb ) to undo the cosf2(r 4 + t [ b ) or sinfi(r 4 + r[ B ) factors before displacement, 
namely, p out = D(f3)R(tt(T^ + r[ B ))p B , we will obtain the best averaged FiQT F+, in the ultraweak coupling limit. 

Mathematically, this can be done by transforming (K B ,A B ) to ( CqK b + tt~ 1 SnA B , CqA b — HSnK B ) in Eq. 
|9| for p B , where Cq = cosf2(r 4 + t[ b ) and Sn = sinfl(r| 4 + t[ b ) [Till . Since the detectors B and C are not 
directly correlated in p AB c> the operation of this counter-rotation on detector B commutes with the joint projective 
measurement on A and C. 

Physically, this may be realized by having Alice continuously send classical signals periodic in her proper time to 
Bob during the whole history, analogous to the local oscillators in optics, so that Bob can determine what rf 4 was 
when the joint measurement on A and C was done, accordingly Bob can counter-rotate detector B for a proper angle 
f^r-j 4 + t[ b ) mod 27 T with t[ b input from his own clock. 


Our numerical results show that this improved protocol is almost the optimal according to (14), though in some 


cases we have to introduce a further squeezing to the coherent state to be teleported in order to optimize the fidelity 
[see Fig. 11 (lower-left)]. 

After introducing the notations and formalism for QT in relativistic consideration, we will now examine carefully 
the special-relativistic effects and the Unruh effect in each of the following four cases. 


IV. CASE 1—ALICE AND BOB BOTH AT REST: TWO INERTIAL DETECTORS 

Let us apply our formulation to the first case, with both Alice and Bob at rest in the Minkowski space and separated 
at a distance d, as the setup in Fig. [T] 


A. Late-time behavior 


The late-time steady state of detectors A and B is simple, in the sense that there is no natural oscillation in time. 
The late-time two-point correlators on the same Minkowski time slice for two UD detectors at rest have been given 
in Eqs. (48)-(51) of Ref. [T 8 |. In these expressions the mutual influences of detectors A and B to all orders (more on 
the mutual influences; see Sec. IV B[ ) are included. From the discussion above Eq. (58) in Ref. JH], one sees that if 
detectors A and B are close enough (d < d en t with the entanglement distance d e nt defined in Ref. pH]), at late times 
these two detectors will have 


((SQ a - SQ b ) 2 )((6P a + 6P b ) 2 ) < U 2 , (33) 

with the operators Q A (t ), P A (t ), Qs(t), and Pg(f) at the same Minkowski time t |43j . This implies that the AB pair 
is in a steady two-mode squeezed state with a phase of 7 r /4 in the Q a Qb subspace of the phase space, and so we may 
be allowed to apply the protocol in Sec. |III| to obtain an averaged FiQT of a coherent state from Alice to Bob or from 
Bob to Alice, 


P _I_i 

1 av ~ j - ^ L Cl — rj 

l + 2h- 1 ^{(6Q.y)((dP + r)/4 


( 34 ) 


in the weak coupling limit according to (14) and beat the classical fidelity F c i. 
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(y = 0.001,0 = 2.3, = 100, d= 1/6} 

{y = 0.001, Q = 2.3, en„«= 100} E N 




T/d 


FIG. 2: (Left) The late-time logarithmic negativity Em (scaled by 3) of two inertial HOs separated at a distance d, with values 
taken on the future or past light cones for one of the two HOs at t = ti I /7 (EnLC, gray solid) and on the ti slice in the 
Minkowski coordinates (EnSM, gray dashed), and the averaged fidelity F av of QT in both teleporting directions subtracted 
by F c i (black). Here, cu max = 100 is the UV cutoff in Eq. ( |35|) . (Right) Em of EnLC evaluated using the cross-correlators 
(7 ?-a(Ii ),Fsiti FT)), 1Z = Q,P at fixed separation d = 1/6 « 0.167, with other parameters the same. (The dotted curve 
represents those negative symplectic eigenvalues that do not count in the definition of Em-) While the two detectors have been 
disentangled according to Em evaluated on the Minkowski time slices (T = 0, EnSM) at this distance, they are still entangled 
around the future and past light cones ( T/d = +1 and —1, respectively; cf. Fig. [I]). 


To look at this possibility more closely, one needs the correlators around the light cone instead of the equal-time 
correlators in the Minkowski coordinates given in Ref. [15] . First, generalize the expressions (52) in Ref. [TBJ to 


^(d,T) S F 


du 


w c cos ujT 


w 2 + 2i^u — LI 2 ± ?F-e iuid 


(35) 


For a given UV cutoff w max , the late-time correlators with detectors A and B at different times, (5Q\(t)) | 7 t»i = 
{SQ 2 B (t + T))| 7t>1 = 2Re[7b+(d,0) + Jb_(d,0)], (SQ A (t)5Q B {t + T ))| 7t »i = 2Re[Jo+(d,T) - Jo_(d,T)], 
(8P%{t ))| 7t »i = (AP 2 (f+ T))| 7t>>1 = 2Re[J- 2+ (d,0) + J- 2 _(d,0)], and (SP A (t)SP B (t + T))| 7t>1 = 2R e[T 2+ {d,T) - 
J- 2 -{d,T)], can be calculated numerically. Using them, one obtains the logarithmic negativity for the EnLC and the 
averaged FiQT between the two detectors by setting t = t-\ and T = ±(d — e) in the above expressions such that 
(t a , t [ b ) = (ti, ty + d — e) or (t[ a , t b ) = (ti, t\ — d + e) and then taking the limit e —► 0+. An example is shown in 
Fig- H It turns out that the late-time EnLC of the AB pair is stronger than the entanglement evaluated on the same 
Minkowski time slice (t a = t b = t\). This implies that the entanglement distance d ent for the EnLC is greater than 
the one we expected according to our old results of entanglement evaluated on the hypersurfaces of simultaneity in the 
Minkowski coordinates (call this EnSM) in Ref. [18]. As one can see in Fig. [2] for the detectors separated at a distance 
d in the range between the entanglement distances for the EnLC and EnSM [0.153 < d < 0.176 in Fig. [2] (left)], the 
averaged FiQT can beat the classical fidelity at late times (F av — F c i > 0) while the detectors are disentangled in view 


of the EnSM (Em = 0 for T = 0). In this range the inequality (141 appears to be violated in view of the EnSM but it 
still holds in terms of the EnLC. Together with the fact that the degree of the EnLC is independent of the choice of 
reference frames and invariant under coordinate transformation, we conclude that the EnLC, rather than the EnSM, 
is essential in relativistic open systems with the “system” consisting of spatially localized objects. 


B. Early-time behavior 

At early times, once Bob enters the future light cone of the spacetime event where detector A started to couple 
to the field, detector B will be affected by the retarded field of A. We call this mutual influence of the first order. 
Detector B will respond to this influence with its backreaction to the field which in turn affects detector A, which 
is called mutual influence of the second order. The subsequent backreaction from A propagates and affects B again, 
which constitutes a mutual influence of the third order, and so on. When the detector-field coupling is not weak 
enough or the spatial separation between the two detectors is not large enough, the higher-order mutual influences 
can get complicated and become very important soon. Fortunately, in the Alice-Rob problem and the quantum twin 
problem to be introduced later, we are working in the weak coupling limit and the retarded distance [31] between 
the two detectors will be very large in most of the history, so the mutual influences are not significant there. To 
compare with those results, assuming that the separation d is large enough, the zero-order result without considering 
any mutual influences in this case would have already been a good approximation at early times. 
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FIG. 3: Spatial and temporal dependence of the logarithmic negativity of EnLC and the best averaged FiQT between Alice 
and Bob. The left plot is for comparison with Fig. 1 in Ref. CHI, with the same parameters there. For the middle and the 
right plots, we set 7 = 0.001, f1 = 2.3, Ao = Ai = 20, r 2 = 1.1, and (d,/9) = (e -ri /\/f2, e _ri \/ff) with n = 1.2. 


t 



FIG. 4: Setup for QT from Alice (thick dashed worldline) at rest to Rob (thick solid worldline) accelerated uniformly from 0 to 
t 2 in his proper time then turning to inertial motion. The hypersurface t = x 1 (blue dot-dashed line) will be the event horizon 
of Rob if r 2 —> 00 . 


We have obtained the evolution in t\ of the logarithmic negativity E^r of the EnLC and the best averaged FiQT 
F+ v in the weak coupling limit, as shown in Fig. 11 (blue curves) for later comparison. The evolution curves are 


roughly exponential decays with small oscillations on top of it at a frequency about twice the natural frequency f l of 
detectors in the weak coupling limit. 

Note that the separation d is also the retarded distance for the classical light signals from Alice to Bob or from 
Bob to Alice. In Ref. El, we see the spatial dependence of the entanglement dynamics: The evolution of quantum 
entanglement between the two inertial detectors, and thus the disentanglement times, depend on d. It is therefore 
not surprising that the evolution of the EnLC and the best averaged FiQT would show a similar dependence on d in 
Fig- [1 The main differences from the EnSM results are the following. First, for the same initial state of the AB pair, 
if the separation d is large enough, one expects that the larger d is, the smaller the “initial” (when r-j 4 = t = 0 + e 
and t[ b = d) EnLC due to the longer time of decoherence of detector B before entering the future light cone of Alice 
emitted at t = e, and thus the shorter is the disentanglement time of the EnLC. Second, the disentanglement rate of 
the EnLC is roughly the same for t < d and t > d, while in Ref. )T5j, we see that the degradation rate of the EnSM 
at early times has nontrivial d-dependence when t > d. 


V. CASE 2—THE ALICE-ROB PROBLEM: ONE INERTIAL, ONE UNIFORMLY ACCELERATED 

DETECTOR 

Our second example has a setup slightly modified from the one in the “Alice-Rob problem” ®HS]. It has been 
claimed that the Unruh effect experienced by Rob (Bob) in uniform acceleration would degrade the FiQT in this 
setup [3]. This is the case in the detector models with the durations of Rob’s constant linear acceleration and the 
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duration of the detector-field interaction being the same and finite, while the teleportation is performed in the future 
asymptotic region when the detectors have been decoupled from the environment |12l . In this section we will examine 
how sound this claim is in our model in which the detectors are never decoupled from the fields and the QT process 
is performed in the interaction region. If Rob is uniformly accelerated, however, there will be an event horizon for 
him, beyond which no classical information can reach Rob (see Fig. [4]). To guarantee the signals emitted by Alice at 
all times can reach Rob to complete a QT from Alice to Rob, we still limit our considerations to the finite duration 
of acceleration, thus no event horizon for Rob, which for all practical purposes is a physically reasonable assumption, 
too. 

Let us consider the setup with Alice at rest along the worldline (f, a -1 — d, 0, 0) with the parameters 0 < (a -1 — d) < 
a -1 and Rob being constantly accelerated in a finite duration 0 < r < f 2 then switched to inertial motion (see Fig. 
0. In the acceleration phase Rob is going along the worldline = (a -1 sinhar, a -1 coshar, 0,0) the same as the one 
for a uniformly accelerated detector with proper acceleration a, and after the moment r = f 2 , or t 2 = a _1 sinhaf 2 
in the Minkowski time, Rob moves with constant velocity along the worldline ((r — r 2 ) coshaT 2 + a -1 sinhaf 2 , (r — 

f 2 )sinhaf 2 + a^ 1 coshaf 2 , 0, 0) in the Minkowski coordinates. Here, the Minkowski time t and the parameter r are 

the proper times of Alice and Rob, namely, t a = t and t b = t. 

Suppose detector C is moving with Alice and its quantum state to be teleported is created right before t = t\, when 
Alice performs a joint measurement on detectors A and C. Then, Alice sends out the outcome carried by a classical 
light signal right after t\ , and Rob will receive the signal at his proper time: 

adv = ^advit x _ / -a -1 In a (a -1 -d-ii) if h < (1 - e~ aT ^)/a - d, . . 

1 1 ( (ti — a -1 + d) e aT2 + a -1 + f 2 otherwise. 

Accordingly, Rob performs the local operation at t b = t'" 1v + e with e —> 0+. 

In the opposite direction, one can also consider the case with detector C moving with Rob, who performs a joint 
measurement on B and C at his proper time t b = n and sends the outcome to Alice by classical channel immediately. 
Then, Alice will receive the message at her proper time, 


.adv _ f adv, \ f d + a 1 (e aTl - 1) if n < t 2 , 

1 [ 11 \d + a ' 1 (e Q A - 1) + (n - f 2 )e Q A otherwise. 


(37) 


and perform the local operation at t a = t\ dv + e. Similar to r adv , here t adv is the advanced time defined by 
| z^(t adv (T)) - 4(r)| 2 = 0 with z° A {t adv {T)) > z%{t). 


A. Dynamics of correlators 


Since Rob stops accelerating at the moment r 2 , the acceleration of detector B is not really uniform. The dynamics of 
the correlators © 0 for nonuniformly accelerated detectors in similar worldlines have been studied in Refs. ESI 132- 
in the weak coupling limit with a not-too-short duration of nearly constant acceleration the behavior of such a detector 
is similar to a harmonic oscillator in contact with a heat bath at a time-varying “temperature” corresponding to the 
proper acceleration of the detector. Analogous to the results in Ref. [32, the dynamics of entanglement here will 
be dominated by the zeroth-order results of the “a-parts” of the self and cross-correlators [3U 131] and the “v-parts” 
of the self-correlators of detectors A and B. The “v-parts” of the cross-correlators are negligible. The higher-order 
corrections by the mutual influences are also negligible in the weak coupling limit with large initial entanglement and 
large spatial separation between the detectors. 

For larger initial accelerations of detector B , the changes of the v-parts of its self-correlators during and after the 
transition of the proper acceleration of detector B from a to 0 are more significant. Consider the cases with the 
changing rate of the proper acceleration of detector B from a finite a to 0 is fast enough so that we can approximate 
the proper acceleration of detector B as a step function of time, but not too fast to produce significant nonadiabatic 
oscillation on top of the smooth variation. According to the results in Refs. m and [31] , for r 2 sufficiently large, the 
v-part of the self-correlators of detector B behave like 


<Wb(t)) 2 ) v 


(Wb(t)) 2 }W+0(t-t 2 ) 


7 ha?e~ 2 ^ T - f J 

67rmo(7 2 + H 2 ) 2 


(((dQ B (oo)) 2 )<°> - ((6Q b ( oo)) 2 )W) 

((SP b (t))\ « ((<5P B (r)) 2 ){“} + 9{t - t 2 ) x 

'{((SP B ( oo)) 2 )W - ((APbM) 2 )i“>) (l - 


(38) 

(39) 
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FIG. 5: Comparison of the best averaged FiQT F+, — F c i (black curves) and the logarithmic negativities Ey (gray) of the 
EnLC from Alice to Rob ](AB), upper plots] and from Rob to Alice [(BA), lower plots], as functions of the moment of the 
joint measurement ti by Alice (ti by Rob) with r 2 = 2 (left) and 10 (right), in the weak coupling limit. Here, a — 1/4 (dotted 
curves), 1/2 (dashed), and 1 (long-dashed gray and solid black). Other parameters are d = 1/4, 7 = 0.0001, 12 = 2.3, Ti = 1, 
77 = 1.2, r 2 = 1.1, ( a,j3 ) = (e~ ri /\/Ti, e~ ri y/Tl), and Ao = Ai = 20. In the upper-right plot, Rob is in the acceleration phase 
when receiving Alice’s signal emitted at ti < (1 — e~ aT2 )/a — d ta 3.42, 1.74, 0.75 for a = 1/4, 1/2, 1, respectively from Eq. 

(36 1 - 


where the superscripts {a} and {0} denote the self-correlators of a UD detector with the same parameters and 
initial state except that it is uniformly accelerated with = a 2 and 0 , respectively, and ((<5Q_B(oo)) 2 )v a ^°^ 


and ((SP B ( oo)) 2 ) 


2\W,{0} 


are those self-correlators in steady state at late times (see Ref. j35|). These approximated 


behaviors have been verified by numerical calculations (see Figs. 3(right) and 4(right) in Re f. |36] ). Note that the 
7 ha 2 term in Eq. (38 1 is actually 0 ( 7 / 0 ), so (( 5Pb(t b )) 2 ) v « Q 2 ((5Qb{ti)) 2 ) v +v, and Eqs. |2lj)-(|24| are still good 
approximations up to 0 ( 7 / 0 ), and we can keep using Eq. (30) here for ro = 0. Below, we apply these approximations 
to calculate the averaged FiQT in the ultraweak coupling limit. 


B. Averaged FiQT in ultraweak coupling limit 


Inserting (ti,r““ w ) in Eq. (36) and (t“ , ti) in Eq. (37) into (t^,t[ b ) in Eq. (25) and its counterpart for the 


opposite teleporting direction, respectively, with the v-parts of the self-correlators (38 1 , (39) and other correlators in 


the approximated form given by Eqs. (21) -(24), we obtain the EnLC and the best averaged FiQT from Alice to Rob 
(E^ B ^ and Fav 3 ^, upper row) and from Rob to Alice ( E and Fa BA ' >+ , lower row) in the sender’s clock in Fig. 
[ 5 ] and in the receiver’s point of view (observed along the past light cones) in the left plots of Fig. [6j respectively. 

The quantities in each plot of Fig. [5]do degrade faster as Rob’s proper acceleration a gets larger and the correspond¬ 
ing Unruh temperature gets higher. However, one has to be cautious at such small accelerations (a = 1/4 to 1 here); 
none of these results can be taken as evidence of the Unruh effect. This is not only because Rob does not accelerate in 
a good part of the histories shown in Fig. [5] but because, more importantly, after the curves in the right plots of Fig. 
[5] are translated to the receiver’s point of view, shown in the left plots in Fig. [6j a larger proper acceleration of Rob 
turns out to give slower degradations of the best averaged FiQT and the EnLC in both teleporting directions even in 
Rob’s acceleration phase. In fact, one can remove the Unruh effect in the calculation by replacing the self-correlators 
of detector B with the Unruh temperature by those for a detector at rest in the Minkowski vacuum, and one will still 
obtain similar curves and the same tendency of the degradation rates against the proper acceleration as those in Fig. 
[5] and the corresponding curves in the left plots of Fig. [6] 
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{r 2 = 10} { £2=2.3, <* = 9.9999} 

E^VlO.2 ^ S) 






FIG. 6: (Left) The black and gray curves are the same results as those in the right plots of Fig. [5] but now against the moments 
Ti dv and ti dv at which Rob and Alice receive the classical signal, respectively. The green and light-green curves represent 
F+ v and Em, respectively, for a = 15 and d = [(2a) 4 + 4 4 ] -1 ' 4 « 0.033. In the up per-left plot when r 2 dv gets large enough 
the curves for the same quantity may cross each other (not shown). From Eq. (371, Alice will receive the signal at t'\ with 
d < t 2 dv < t adv (f 2 ) « 44.98, 295.08, 22025.7 for a = 1/4, 1/2, and 1 if Rob emits the classical signal in his acceleration phase. 
(Right) E^ b) and E^ A) for the EnLC at fixed moments r“ d ” = 9.9999 and ti dv = 200 in Rob’s and Alice’s points of view, 
respectively, as functions of a (black). The gray dotted curves are the same quantities with the Unruh effect removed from the 
self-correlators of detector B. Here, d = [(2a) 4 4- 4 4 ]^ 4 ' 4 , and t 2 = 10, so that, in the lower-right plot, if a > 0.45, Rob will be 
in the acceleration phase when he performs the joint measurement as the sender. Other parameters are the same as those in 
the previous figure. 


The behavior of the curves in Fig. [5] can be explained simply by the go-away setup in the Alice-Rob problem 
and the Doppler shift. For Fa AB ^ + and from Alice to Rob with t\ and f 2 fixed, the proper time T ddv in Eq. 

(36) when Rob receives Alice’s signal increases rapidly as the value of a increases, which allows for a much longer 
duration of decoherence for detector B before Rob’s operation. This yields a higher degradation rate in t\ (Alice’s 
clock) for larger a in the evolution of the best averaged FiQT from Alice to Bob. On the other hand, Alice’s signal 
is more redshifted and so Alice’s clock looks slower for a larger a in Rob’s point of view. When a is not too large, 
the apparent slowdown of decoherence for detector A can beat the increasing rate of decoherence time for detector B 
such that the larger a is, the slower is the degradation in T° dv [see the black and gray curves in Fig. 6] (upper-left)]. 
Similarly, for a fixed value of a, Eq. (|36|) implies that T ddv for Rob grows rapidly as the duration of Rob’s acceleration 


phase r 2 increases, which causes a much faster degradation of F^ 3 ^ and E^ 3 ' 1 in t\ also. Indeed, the curves in the 
upper-right plot (f 2 = 10) of Fig. [ 5 ] drop faster than those in the upper-left plot (f 2 = 2) for each value of a. Let t c i 
be the moment of t\ when Fay 3 ^ drops to the value F c i for the classical teleportation. When af 2 is sufficiently large, 
T adu k e g0 i ar g e that f cl a~ l — d, which is the moment in Alice’s clock when Alice crosses the event horizon for 

Rob as r 2 —> 00 . For Fa BA ^ + and E in the opposite teleporting direction, the situations are similar, even though 
ostensibly there is no event horizon for Alice. 

This is not the whole story, though. If we increase Rob’s proper acceleration a further, while the EnLC from Rob 
to Alice Ecj BA ^ is always an increasing function of a [Fig. [i] (lower-right)], such a tendency will be altered when 
a > O(fl) in the EnLC from Alice to Rob E^ B \ as shown in Fig. [^(upper-right), mainly by the factor coth(7rQ/a) 
in the self-correlators of detector B , e.g., ((^(r)) 2 )^ « (h/2tt) coth(7rfl/a)(l — e~ 1T ), for Eq. (22) when Bob is 
accelerated m ■ Only in this regime, the Unruh effect is significant and dominates over the apparent slowdown of 
Alice’s clock observed in Rob’s acceleration phase, in the sense that a higher Unruh temperature leads to a higher 
degradation rate of the best averaged FiQT and the EnLC. After Rob’s acceleration phase is over, however, due to 
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the higher relative speed between Alice and Rob causing a stronger redshift of Alice’s clock signal with a larger a, the 
degradation later in Rob’s point of view can be slower than those with a smaller a. Indeed, we see that the slopes 
of the black and gray dotted curves (a = 1/4) are more negative than the slopes of the green and light-green curves 
(a = 15), respectively, for T° dv > f 2 = 10 in Fig. [ 6 ] (left). 

Comparing the upper and lower plots in Fig. T5j one sees that the moment t c i (or t c i defined similarly for Rob) 
when QT from Alice to Rob (or from Rob to Alice) loses advantage over “classical” teleportation is always earlier 
than the disentanglement time evaluated around the future light cone of the joint measurement by Alice (or Rob) at 
t-[ (or Ti). This confirms that the EnLC of the AB pair is a necessary condition for the best averaged FiQT beating 
the classical one, as indicated in Eq. (14). 


C. Beyond ultraweak coupling limit 

Beyond the ultraweak coupling limit, both the averaged fidelity F av and the logarithmic negativity Ey are strongly 
affected by the environment. In the cases in which mutual influences to the first few orders are small compared with 
the zeroth-order, quantum entanglement of detectors A and B disappears quickly due to the strong corrosive effects 
of the environment. We expect that the best averaged fidelity F^ 3 ^ and Fa 3A ^ + would drop below F c i even quicker 
[HI- Similar results on entanglement were given earlier in Ref. |19j , though the degrees of entanglement in Ref. [19] 
are evaluated on the time slices in the Minkowski coordinates or Rindler frames rather than those evaluated around 
the light cones. 


VI. CASE 3—QUANTUM TWIN PROBLEM 


In the above results, we have seen that the relativistic effects entering the description of the dynamics of the 
detector pair can dominate over the Unruh effect experienced by the accelerated detector in the degradation of the 
best averaged FiQT and the EnLC between the pair. The apparent “slowdown” in the dynamics of the sender in the 
viewpoint of the receiver in a QT process can be perceived by the receiver in the redshift of the clock signal from the 
sender. Nevertheless, in the setup of the Alice-Rob problem, since the retarded distances from Alice to Rob and from 
Rob to Alice are always increasing in time, only the redshift of the clock signal from the other would be observed, and 
so both Rob and Alice would conceive that their partner’s clocks are always slower than their own. One may wonder 
what will happen when Alice and Rob (Bob) undergo more general motions. 

To get a more comprehensive picture, a simple but helpful extension is to consider a setup similar to the classical 
twin “paradox” Wl in which we would have a consistent description of the asymmetric aging, red- and blueshifts of 
the clock signals, and the inertial and noninertial motions. Indeed, recall that in special relativity the twin paradox 
originates from the disparity between Alice the twin at rest and Bob the traveling twin: Alice seeing Bob going away 
is the same as Bob seeing Alice going away, so each one is supposed to observe the other with the same time dilation. 
Why does Bob become younger but not Alice when they meet again? The resolution is that, for Bob to return to 
Alice, he must turn around at some point, thus undergoing some period of acceleration, and the principles of special 
relativity do not apply to noninertial frames. When coupled to quantum fields, the Unruh effect experienced by Bob 
during the periods of acceleration will come into play. With the theoretical tools developed and knowledge gained in 
the previous sections, luckily, this quantum twin problem becomes straightforward. 

Suppose Alice is at rest with the worldline 4 = ( t , — d, 0,0), d > 0 and the proper time t a = t, Bob is going along 
the worldline z b (t) with z B = z B = 0 and 


(4(4 4(a)) 


' (u °) 

(i sinha(r - fi) +fi, ^[cosha(r - f 4 ) - 1]) 

(72(T-f 2 )-|-4(f 2 ), 72^2(T-f2) + 4(f 2 )) 

< (i[sinha(r - r 3p ) - 72 i; 2 ] + 4 (t 3 ), 4 [cosh a(r - r 3p ) - q 2 ] + 4 (t 3 )) 
( 72 ( 1 - - f 4 ) + 4(f 4 ), - 72 ^ 2 (r -t 4 ) + 4(7 4 )) 

(a [sinha(r - r 5p ) - -/ 2 v 2 } + 4( r s)> i[cosha(r - r 5p ) - y 2 ] + 4( T 5)) 

. {(t - T 6 ) + z%(f . 6 ), 0 ) 


0 < r < t 1; 
fj < r < t 2 , 
t 2 <t < f 3 , 
for f 3 < r < f 4 , 
f 4 <T < f 5 , 
T 5 < T < f 6 , 

T > 4 


(40) 


where t b = t is Bob’s proper time, r p = r 2 - r 4 = (r 4 - t 3 )/2 = t 6 - r 5 , r 3p = t 3 + r p , r 5p = t 5 + r p , t 3 - r 2 = r 5 - r 4 , 
72 = coshaf p , and r y 2 v 2 = sinhaf p (see Fig. [T| . Here we set the minimal distance between Alice and Bob d to 
be sufficiently large to avoid the singular behavior of the retarded fields, and thus the mutual influences, when the 
detectors are too close to each other in the final stage (for example, see Ref. jT5] ). 
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FIG. 7: QT between Alice (thick dotted worldline) and Rob (thick solid worldline) in a setup of the twin problem, where the 
worldline of the traveling twin Bob is given in Eq. (|40[). 


A. Evolution of correlators 

Below, we consider a case in the ultraweak coupling limit, with Bob still at his youth (yfg -C 1) at the moment 
when he rejoins Alice, who is also in her early age ( 72 ° (fc) < 1) but much advanced in age than Bob at that moment 
[e.g., f 6 = 16 for Rob and z°(t 6 ) = 220 for Alice in Figs. |8]and[9]. 

As before, suppose the combined system is initially in a product state ®Pab ®Pc to P °f th e well-studied 

self-correlators for a detector at rest in Minkowski vacuum [35] , the subtracted v-parts of the self-correlators of detector 
B USUI] in our weak coupling limit, S(TZ b (t), 1Z b (t)) v = {'R- b (t),1Z' b (t)) v — (7^ s (r), 7^^(r)) v | a=0 , TZ,TZ' = SQ,SP 
have been obtained numerically. We found that 5(1Z, TZ’) V starts to oscillate after the launch of Bob. The oscillations 
would be amplified whenever the acceleration suddenly changes from one stage to the next due to the nonadiabatic 
effect m, while its mean value grows due to the Unruh effect when detector B is undergoing accelerations and decays 
during the time intervals in the inertial motion. Anyway, the amplitude of 5(TZ B ,1Z B ) V is always as small as O(y) 
compared with (1Z B ,1Z' B ) V1 while ( 1Z B ,1Z' B ) V is small compared with (lZ B ,lZ' B ) a in such an early stage. 

We further obtained the numerical results for the cross-correlators between A and B, (TZ A (t), TZ' B {r adv {t) — e)), 
and ( lZ A {t adv {r ) — e),TZ' B (T)), around the future light cone of Alice and Bob at t a = t and t b = r, respectively. 
We find that they oscillate in time about zero during the whole journey of Bob until he meets Alice again. The 
oscillations appear irregular since the motions and the time dilations of the two detectors are asymmetric. While 
the amplitudes of the oscillations of the a-parts of the cross-correlators are 0(1), the amplitudes of the v-parts are 
0 ( 7 ) and negligible in the weak coupling limit. After Bob returns and both detectors are at rest, the behavior of the 
a-parts of the cross-correlators continues in the same way, but the v-parts of (Qa,Qb) and ( Pa,Pb) turn into small 
oscillations on top of slow growths or decays in time, similar to those in the cases with two detectors at rest (see Sec. 
V in Ref. [IB]). 

Corrections from the mutual influences (7£- 0 \ i,j = A, B up to the first order of 7 /d have been worked 

out to check the consistency of our approximation. There is one correction to each of both the a-part and v-part 
of the correlators (Qf) and (Pf) and two corrections to those for the other correlators. Thus, we have a total of 32 
corrections of the first order. We find that, during Bob’s journey, the corrections to the v-part and the a-part of each 
correlator are 0 ( 7 ), small compared with the zeroth-order results. After Bob returns and stays at rest by Alice, these 
corrections from the mutual influences start to grow in magnitude. If the separation d of Bob and Alice is small, 
these corrections may overtake the zeroth-order results and one has to include higher-order mutual influences |18j . 
Here, we simply terminate our simulation at t b = ff ss 24 in Bob’s proper time, which is early enough to justify our 
first-order approximation. 

One may worry that the backreaction from detector B to the field during r £ (r 4 , r 6 ) would form a shock wave and 
hit detector A in the period when Bob heads back to Earth and decelerates [t £ {t a dv{f 4 ),tadv(f 6 )) ~ (220.88,221.43) 
in the left plots of Figs. [ 8 ] and [9], analogous to the shock electromagnetic wave along the past horizon of a uniformly 
accelerated charge in classical electrodynamics [35]. Fortunately, in our results, these mutual influences do not 
significantly impact on detector A since they are off resonant. 


B. Entanglement dynamics 


With the results of the correlators we are able to calculate the dynamics of the EnLC in both teleporting directions. 
Our first example is shown in Fig. [s] In the left plots, one can see similar decays of E^ B ' 1 (corresponding to the QT 
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FIG. 8: Dynamics of the EnLC in both teleporting directions including first-order correction from the mutual influences, in the 
clocks and points of view of Alice (upper-left and lower-left plots) and Bob (lower-middle and upper-middle plots). The gray and 
pink regions in these plots represent the three time intervals when Alice’s signal reaches Bob or Bob’s signal is sent to Alice during 
Bob’s acceleration phase (the leftmost diagrams). Here, 7 = 0.001, If = 2.3, and Ao = Ai = 20, a = 2, (a, (3) = (1.4, 0.2), 
and the initial or final spatial separation d = 1. For Bob’s worldline, we set (to, 77 , f 2 , f 3 , f 4 , f 5 , f 6 ) = (0,2,4,7,11,14,16). 
Other parameters have the same values as those in case 2. (Upper right) AE^ B ' 1 = E^ B ' 1 — E^ BA0) is the correction of 
entanglement strength to the zeroth-order result from the mutual influences up to the first order. (Lower right) 

A E (BA) = E^ A) - E^ A)(0) is similar. 


from Alice to Bob) in Alice’s clock and e\^ A ^ (from Bob to Alice) in Alice’s point of view. While in the middle plots 
the two curves in Bob’s clock or point of view drop significantly in different periods, the values of E^ B ' > and E^ A ^ 
around the moment when Bob comes back to Alice are quite the same. Once again, the details of the history depend 
on the point of view, but here we further see that different views on the EnLC tend to agree when Bob rejoins Alice. 
The reason is simple. When two detectors are close enough, the amplitudes of the mode functions in the operators 
Q, P of detectors A and B at t a = t and t b = T adv (t), respectively, are relatively close to the ones at t a = t adv {r) 
and t b = r if d <C c/ 7 . So, these operators give similar expectation values of the two-point correlators with respect 
to the same initial state. In the case Rob never returns, as in the Alice-Rob problem studied in the previous section, 
E^ B ' > and E^ A ^ in different teleporting directions will never be commensurate after the initial moment. 

In our example, the mutual influences tend to enhance the entanglement during the space journey of Bob. Denote 
the zeroth-order results of the logarithmic negativities for the EnLC as E^-\ and the enhancement by the mutual 
influences as A Ejj- = (Ey — Ey' 1 ). In the right plots of Fig. [ij we find both A Ey 8 ' 1 and A Ey A ^ grow from zero to 
some value when Bob launches (r, r adv ss 77 ), and then during Bob’s journey A E ( y B ' 1 and A Ey A ^ roughly remain 
constant between +0.0014 to +0.002, which is of the same order as 7 /d ~ 10” 3 . However, when Bob returns to Alice, 
the corrections to the logarithmic negativity from the mutual influences oscillate between positive and negative values 
with the amplitudes increasing in time. 

Furthermore, in the right plots of Fig. [ij one can see that A Ey' B ' 1 appears to be slightly “kicked” at about 
t G {t a dv(f 4 ),tadv(T 6 )) « (220.88,221.43) and A E^ A ^ at about r « 15 G (t 5 ,t 6 ). This could be due to the shock 
waves emitted by detector B during r G (t 4 ,t 6 ) that all hit detector A at t « 221. In our results the impact of the 
first-order correction never gets significant compared to the zero-order correlators. 

C. Quantum teleportation 

Next, to compare the averaged FiQT, we set ( a,j3) = (e -ri y/h/fi, e~ ri \fhVL) , = 1.2 for the initial state of the 

entangled pair of the detectors as the one in the previous section. The results are shown in Fig. [9] Again one can see 
that the evolutions of the best averaged FiQT F+, in either teleporting direction subtracted by F c i is similar to the 
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FIG. 9: The averaged FiQT of a coherent state of detector C from Alice to Bob (Fa AB ^) and from Bob to Alice with 

(black curves) and without (purple) using the improved protocol in the viewpoints of Bob (left) and Alice (right), respectively. 
Here the entangled pair starts initially with (a, ft) = (e~ ri /VTl, e~ ri \/9), fl = 1.2, and we assume the joint measurements of 
detectors C and A by Alice or C and B by Bob collapse the measured detector pair to a squeezed state with squeeze paramater 
T 2 = 1.1. Other parameter values are the same as in the previous figures. The scaled logarithmic negativities of EnLC with 
the same parameters are plotted in blue curves for comparison. One can see that the evolution of Ej y in time is similar to 
Fav ~ Fd- 


evolution of the logarithmic negativity Eff of the EnLC of detectors A and B. 

We keep the curves for the averaged fidelities F av without using the improved protocol in the upper row of Fig. [9] to 
give the readers a flavor how the sender’s clock is observed by the receiver [recall Eqs. (301 and (31)]. One can see that 
there is no significant enhancement of decay for F+ v or Ey due to the Unruh effect when Bob is in any acceleration 
phase (gray or pink regions), since we take the proper acceleration a = 2 for Bob, which is not too large there. In 
contrast, significant drops of or E^y in Fig. [9] (left) happen between the second and the third acceleration 

phases, when Bob sees a strong blueshift in the clock signal emitted by Alice and so Alice’s clock looks much faster 
than Bob’s in his viewpoint during this period [Alice’s signal emitted during (t 4 ,t 5 ) = (28.836,192.63) reaches Bob 
during the period (r adu (f 4 ), r adv (t 5 )) = (f 4 ,r 5 ) = (11,14)]. This implies that quantum coherence of detector A in 
this period fades much more quickly than any other period in Bob’s viewpoint so that quantum entanglement and 
the best averaged FiQT are degraded faster in this stage. The significant drops of the EnLC in the middle plots 
of Fig. [s] are due to the same reasons. For Fay A ^ + and E in Fig. [9] (right), the drop is much less significant, 
though. This is because the period in which Alice receives similar blueshift clock signal from Bob is much shorter 
than the time scales of decoherence (I /7 = 1000) either in Bob’s clock (r 5 — r 4 = 3) or in Alice’s point of view 
(t adv (f 5 ) - t adv (t 4 ) = (f 5 - f 4 )e-°^ « 0.055). 

In the above cases we have seen that the relativistic effects play a dominant role in QT. One can ask when the 
Unruh effect will become more significant in the QT from Alice to Bob. Our results so far show that this happens 
only in Bob’s point of view and only when Bob’s proper acceleration a is large enough (see Figure [6] (upper row), for 
example). In other words, only in a highly accelerated receiver’s point of view can this happen. One can construct 
setups in which the Unruh effect can be singled out, such as those with both detectors uniformly accelerated or in 
alternating uniform acceleration [Fig. 10 (middle and right)], but then the receiver is also accelerated in these setups 
after all. Is it possible for a receiver in QT remaining at rest to see the domination of the Unruh effect? With this 
aim, we construct below a setup with Alice at rest, while the relativitic effects of time dilation and varying retarded 
distance are suppressed and the Unruh-like effect are significant in QT in both directions. 


VII. CASE 4—TRAVELING TWIN IN ALTERNATING UNIFORM ACCELERATION 


To highlight the regimes in which the Unruh effect stands out in comparison with other relativistic effects, we design 
a case in which Bob the traveling twin undergoes an alternating uniform acceleration (AUA) considered in Ref. |26| 
with the period of motion so short that the maximum speed of Bob is low enough and the retarded distance between 
Alice and Bob does not vary too much, while the proper acceleration can still be very high. Consider the case with 
Alice at rest along the worldline (t, —d, 0,0) and Bob going along the worldline 


(''■)= - 


sinho ( r — n— 
2 


+ 2 n sinho— 
4 


(-ir 


cosh a [ r — n-y ] — cosh a-y 


, 0,0 


(41) 


with n(r) = Floor} (2r/r) + (1/2)}, linearly oscillating in the x 1 axis about the spatial origin [see Fig. 10 (Left)], 


where r p is the period of Bob’s oscillatory motion in his proper time. In this case the classical light signal emitted by 
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FIG. 10: (Left) QT from Alice (thick dotted) to Rob (thick solid) in a variation of the twin problem, where the traveling twin 
Bob is in alternating uniform acceleration with the worldline (411. One can conjure up settings that single out the Unruh effect, 
such as letting both Alice and Bob be uniformly accelerating (middle) or both in alternating uniform acceleration (right), where 
n is an integer. Note that in the middle plot the relativistic effects in affecting the description of the dynamics are totally 
suppressed only in the one-way QT from Alice to Bob, but not from Bob to Alice. 


Alice at t will reach Bob at 


r adv {t) = n-E - 

w 2 a 


log |cosha-^ + (— l) r 


2nsinha— — a(t + d) 
4 


(42) 


where h(t) = Floor}(2t/f p ) + (1/2)} with t p = 4a 1 sinh(ar p /4), while the classical light signal emitted by Bob at 
r will reach Alice at t adv (T) = d + z° b {t) + z\ (r). To compare with cases 2 and 3 in which the mutual influences 
are small, the retarded distance between Alice and Bob is set to be large enough. Also when the period of motion 
is much less than the natural period of the detector (f p <C T = 2n/fl), the time-averaged subtracted Wightman 
function will be a good approximation in calculating the self-correlators of detector B (see Sec. 5.1 in Ref. ffi?]). 
These assumptions simplify the calculat ion very much in the weak coupling limit. 

We show some selected results in Fig. |TTJ For the logarithmic negativity E^ B ' 1 of the EnLC and the best averaged 
FiQT Fav B ^ + from Alice to Bob in Alice’s clock or in Bob’s point of view, when a is small and f p is large, the 
disentanglement time for the EnLC of the joint measurement by Alice is still longer than the one in case 1 with the 
same parameters except a = 0. Here, time dilation of detector B dominates. As a gets larger, with the maximum 
speed fixed ( af p =constant), one starts to see the evolution curves for E^ B ' > and drop faster than the ones 

with a = 0 in some parameter range of (a, (3) for the initial state ([3| [Fig. |TT] (Left) in Bob’s point of view; the 
plots in Alice’s clock look similar]. When a is large enough, the initial states with all values of (d,/9) will see faster 
degradations of the EnLC and the best averaged FiQT, both in Alice’s clock or in Bob’s point of view than those in 
the a = 0 case [Fig. [IT] (Middle)]. Now, we can say that the Unruh effect dominates, though the effective temperature 
experienced by detector B is lower than the Unruh temperature with the averaged proper acceleration a |26j . In 
the reverse teleporting direction, for the logarithmic negativity E^ A ^ of the EnLC in Alice’s point of view, we see 
clearly that the larger a is, the shorter the disentanglement time in Fig. [ll] (right), where the Unruh effect has been 


dominating the degradation of E from a = 10 for all values of (a,/3), while EYf a> with a = 10 still has a longer 


(AB) 


disentanglement time than the one with a = 0 in a corner of the parameter space around (a, /3) 
in the lower-left plot of Fig. US- 


(1.4, 0.2), as shown 


One interesting observation in calculating Fig. 11 (lower left) is that when a is large enough the averaged FiQT of 


a coherent state using the entangled AB pair initially with (a,/3) in some finite parameter range will never achieve 
Eav or F av A F cl = 1/2. One has to modify the quantum state to be teleported from a coherent state to a 
squeezed coherent state with the squeezed parameter ro > 0 in Eq. ([2]) and tune the value of ro to push the averaged 
FiQT above F c j toward the optimal fidelity F opt in (141, so that the time t c i when F av — F c i touches zero is closer to 
the disentanglement time tdE of the EnLC. Note that ro itself is a-part of the protocol and not among the quantum 
information to be teleported. 


VIII. SUMMARY AND DISCUSSION 

We have considered the quantum teleportation of continuous variables applied to three Unruh- DeWitt detectors 
with internal harmonic oscillators coupled to a common quantum field. The basic properties of relativistic effects in 
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FIG. 11: Dynamics of the EnLC and the FiQT between Alice and Bob with Bob at rest (blue curves), in AUA (gray and 
black), and as in the twin problem (purple). The mutual influences are ignored, and the initial state of the AB pair has 
(a ,/?) = {e~ ri e~ ri \/fi) with n = 1.2 (upper row), or (1.4, 0.2) (lower). (Left) The scaled E^ B ' 1 of the EnLC (lighter) 
and Fav BS>+ (darker) from Alice to Bob subtracted by the classical fidelity F c i = 1/2, with d = 1 both for Bob in AUA and 
at rest, in Bob’s point of view. Bob in AUA has a = 10 and the period of his oscillatory motion f p = T/ 16, T = 2'k/Q,. The 
squeezed parameter in p'^c is r 2 = 5.1, and other parameters are the same as before. In the lower-left plot the teleported state 
has ro = log 2. (Middle) Comparison of the EnLC between Alice at ti and Bob at T adv (ti) in different motions in Alice’s clock. 
Here, d = 4, a = 2 in the twin problem and a = 20 in the AUA case in which t p = T/ 32 for Bob. (Right) Dynamics of the 
EnLC between Bob at n and Alice at ti = t adv (r 1 ) in Alice’s point of view, where Bob is at rest (blue, a = 0) or undergoes 
AUA (gray, from dark to light a = 2 n ■ 10, n = 0 to 7 with af p = 10T/16 fixed). Again d = 4 with other parameters unchanged. 


dynamical open quantum systems such as the frame dependence of quantum entanglement, wave functional collapse, 
Doppler shift, quantum decoherence, and the Unruh effect have all been considered consistently and their linkage 
manifestly displayed. Below is a summary of what we have learned from these studies. 


A. Entanglement around the light cone 


Quantum entanglement of two localized objects at different positions requires the knowledge of spacelike correlations, 
while the averaged FiQT involves timelike correlations between two causally connected events. In general these two 
quantities are incommensurate. To compare them, in Sec. m we introduced the projection of the wave functional 
around the future light cone of the joint-measurement event by the sender, so that right after the wave functional 
collapse the sender’s classical signal of the outcome reaches the receiver, according to which the receiver performs the 
local operation immediately. The averaged FiQT obtained in this way can be directly compared with the degree of 
quantum entanglement in the entangled detector pair evaluated right before the wave functional collapse, namely, the 
EnLC, which can be easily calculated in the Heisenberg picture. 

We have observed that the best averaged FiQT always drops below the fidelity of classical teleportation earlier than 
the disentanglement time for the EnLC in each of our numerical results. This confirms the inequality , which 
implies that entanglement of the detector pair is a necessary condition for the averaged FiQT beating the classical 
fidelity. In Sec. IV A we further showed that the inequality (141 may appear to be violated by the degrees of quantum 
entanglement evaluated on a time slice in conventional coordinate systems. This proves that the EnLC, rather than 
the conventional ones, is essential in QT in a relativistic open quantum system. 


B. Multiple clocks and points of view 

For a relativistic system including both the local and nonlocal objects such as a detector-field interacting system, 
the Hamiltonian, quantum states, and quantum entanglement extracted from the states all depend on the choice of 
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the reference frame |19j . Part of the coordinate dependence can be suppressed by evaluating the physical quantities 
around the future or past light cones of a local observer. However, this does not give a unique description on a 
physical process, since each local object has a clock reading its own proper time, which is invariant under coordinate 
transformations. In particular, a QT process involves two different physical clocks for the sender and the receiver 
localized in space, and the degradation of the EnLC and the averaged FiQT in the same process can appear very 
differently in the sender’s clock and in the receiver’s point of view along his/her past light cone. When describing 
nonlocal physical processes with local objects in a relativistic open quantum system, one has to first specify which 
clock or which point of view being used; otherwise, there will be ambiguity in the statements. 


C. time dilation, Doppler shift, and acceleration 

It is easy to understand that the FiQT between localized quantum objects in a field vacuum with one party or both 
accelerated would be degraded by the Unruh effect because of the thermality appearing in these accelerated objects 
[5]. However the more ubiquitous relativistic effects in inertial frames such as time dilation and Doppler shift (related 
to the relative speed) that are mixed in with effects due to acceleration have not been understood fully in the context 
of QT. These effects and their interplay are the focus of this study. What we found that may be surprising is that the 
relativistic effects in affecting the description of the dynamics can overwhelm the Unruh effect. For example, there 
is degradation of fidelity when both parties are inertial, as shown in our case 1, and a larger acceleration does not 
always lead to a faster degradation, as shown in our cases 2 and 3. 

The averaged FiQT in cases 2, 3, and 4 do depend on the proper acceleration a in Bob’s acceleration phase 
significantly. In case 2, we find that the larger a is, the higher the degradation rate will be in the sender’s clock for 
the best averaged fidelities F+, of QT both from Alice to Rob and from Rob to Alice. Nevertheless, the increasing 
redshift as the retarded distance between Alice and Rob increasing indefinitely in time is the key factor for the a 
dependence here. In the receiver’s point of view, that the degradation rate increases as a increases, is true only for a 
receiver accelerated with proper acceleration large enough, when the Unruh effect fully dominates. In case 3 a larger 
a turns out to give a longer disentanglement time of the EnLC in the clock of the sender Alice at rest. The key factor 
there is that detector B with the traveling twin Bob ages much slower than detector A with Alice at rest when they 
compare their clocks at the same place after Bob rejoins Alice. The acceleration of Bob leads to this asymmetry of 
time flows as in the well-known twin paradox and Bob’s slower clock helps to keep the freshness of quantum coherence 
in the AB pair longer from the view of Alice’s clock, while the retarded distance between Alice and Bob is bounded 
from above. 

To suppress the relativistic effects in what is observed by Alice, who is always at rest, we considered case 4 in 
which Bob is undergoing an alternating uniform acceleration with a small speed and a constant averaged retarded 
distance. The results indeed show that the larger the a, the shorter the disentanglement time for EnLC, even in 
Alice’s point of view when a is large enough, although the Unruh temperature is not well defined in this setup for 
the lack of a sufficiently long duration of uniform acceleration. 
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Appendix A: Reduced state of a detector with its entangled partner being measured 


In our linear system the operators of the dynamical variables at some coordinate time x° = T of an observer’s 
frame after the initial moment Tq are linear combinations of the operators defined at the initial moment |34j : 


Q d (r d (T)) = Y, [$ (r d )Q l $ + ./d d V d )P™ 

d' 

Km = Y [** cn^d' + 


H J d 3 y [^(r d )f+/ d y (r d )f 

d 3 y k(T)4M + 


(Al) 

(A2) 


from which the conjugate momenta P d (T) and n x (T) to Q d {T) and $ X (T), respectively, can be derived according 
to the action jlj). Here we denote = 0^(T n ) (e.g., $y' = <l(T ra ,y) and fly ' = n(T„,y)), and all the “mode 
functions” <jr^([T) and /|(T) are real functions of time (£,£, v £ {A,B,C} U {x}, x £ R 3 in (3+l)-dimensional 


Minkowski space), which can be related to those in k space in Ref. [33j. Then from Eqs. (Al) and (A2), those 
correlators in Eqs. 0 0 can be expressed as combinations of the mode functions and the initial data, e.g., 


( Qa( t a )) — < I ) a( t a) ( I ) a{ t a)((Qa ) 2 )o + 

[ d 3 xd 3 y^(r A )^(r A )($[ £ ° 1 ,l>M ) 0 


(A3) 


where (• • -) n denotes that the expe ctatio n va lues are taken from the quantum state right after x° = T n . 

Comparing the expansions (Al) and (A2) of two equivalent continuous evolutions, one from x° = Tq to x° = T) 
then from x° = T\ to a; 0 = T2 and the other from x° = Tq all the way to x° = T2, one can see that the mode functions 
have the identities, 
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0f 1] ^[i°] + f" [21] 4 [10] , 
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where the DeWitt-Einstein notation with v £ {A, B, C} U {x} is understood, F^ mn 1 = F(T m — T n ), and ^ d (T d (T)) = 
( rd ( T ))> ni(T) = d 0 (l)i{T), p C d {r d (T)) = 9 d /^(r d (T)), and p^(T) = ^ 6{T) in the momentum operators. 
Similar identities for 7r^ and p ^ can be derived straightforwardly from Eqs. (A4) and (A5). Such identities can be 
interpreted as embodying the Huygens principle of the mode functions and can be verified by inserting particular 
solutions of the mode functions into the identities. 

In Ref. [23] one of us has explicitly shown that in a Raine-Sciama-Grove detector-field system in (l+l)-dimensional 
Minkowski space, quantum states in different frames, starting with the same initial state defined on the same fiducial 
time slice and then collapsed by the same spatially local measurement on the detector at some moment, evolve to the 
same quantum state on the same final time slice (up to a coordinate transformation), no matter which frame is used 
by the observer or which time slice is the wave functional collapsed on between the initial and the final time slices. 
This implies that the reduced state of detector B at the final time is coordinate independent even in the presence 
of spatially local projective measurements. For the Unruh-DeWitt detector theory in (3+l)-dimensional Minkowski 
space considered here, the argument is similar, as follows. 

Right after the local measurement on detectors A and C at Xj (for a simpler case with the local measurement only 
on detector A, see Ref. PH), the quantum state at Tj collapses to Pac®Pb<s> on the Tj slice of the observer’s frame. 
Similar to Eq. (|9j), here Pb<s>^ f° r detector B and the field $ x in the postmeasurement state is obtained by 
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where p is the quantum state of the combined system evolved from Tq to Tf and cr £ {H}U{x}. Since p AC is Gaussian, 
a straightforward calculation shows that p B $ has the form 
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Again we use the DeWitt-Einstein notation for </,£ £ {B} U {x}, which run over the degrees of freedom of detector 
B and the field defined a t x on the whole time slice, n running from 1 to 4 corresponds to the four-dimensional 


Gaussian integrals in Eq. (A6). depends only on the two-point correlators of detectors A and C at the moment 


of measurement, while (T^) and A4i ra; (II^) are linear combinations of the terms with a cross-correlator between 

detector A or C and the operators or 11^ (4 >b = Q B and n B = P B ), respectively, multiplied by a few correlators 
of A and/or C , all of which are the correlators of the operators evolved from Tq to Tf with respect to the initial state 
given at Tq. This implies that the two-point correlators right after the wave functional collapse become 
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For example, (( SQ [ B ] ) 2 ) 1 = Q BB (T{) - 


(n) 


refers to the operator O b defined at Ti and refers to the operator O b (Ti — T 0 ) in the Heisenberg picture. 

Suppose the future and past light cones of the measurement event by Alice at x° = T± crosses the worldline of Bob 
at his proper times and r[ et , respectively. At some moment in the coordinate time a; 0 = T M of the observer’s 
frame before detector B enters the future light cone of the measurement event on detector A, namely, when Bob’s 
proper time t b = t(T m ) £ ( r [ et , r{ ldv ), the two-point correlators of detector B are either in the original, uncollapsed 
form, e.g., (( SQ B ) 2 (T M — T 0 )) 0 , if the wave functional collapse does not happen yet in some observers’ frames, or in 
the collapsed form evolved from the postmeasurement state, e.g., 
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in other observers’ frames. Here we have used the Huygens principles (A4) and (A5) and defined 
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with &a,c = Qa,c and fl a,c = Pa,c, while I^ is derived from those J^ n> and J c 
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that before detector B enters the light cone, one has <f> 


A[M1] 
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_ ,A[M1] 
— Jb 
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= 0 , 


such that Y^ /0 ^ reduces to Q l L 


pairs in Eqs. (A8)-(A10). Note 

jmf 


So at 


the moment T M . the correlators of detector B do not depend on the data on the T-\ slice except those right at the 
local measurement event on detectors A and C. This means that, once we discover the reduced state of detector B 
has been collapsed, the form of the reduced state of B will be independent of the moment when the collapse occurs 
in the history of detector B (e.g., t b = t b or r'f in Fig. [lj, namely, the moment at which the worldline of detector 
B intersects the time slice that the wave functional collapsed on. 

No matter in which frame the system is observed, the correlators in the reduced state of detector B must have become 
the collapsed ones like Eq. (All) exactly when detector B was entering the future light cone of the measurement 
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event by Alice, namely, t b = , after which the reduced states of detector B observed in different frames became 

consistent. Also after this moment, the retarded mutual influences reached B such that q and /'^ M1 ^ would 
become nonzero and get involved in the correlators of B. In fact, some information of measurement had entered the 
correlators of B via the correlators of A and C at t\ at the position of Alice in and much earlier. 

Nevertheless, just like what we learned in QT, that information was protected by the randomness of measurement 
outcome and could not be recognized by Bob before he has causal contact with A lice. 

Thus, we are allowed to choose a coordinate system with the T M in Eq. \A1\ I giving t b (T m ) = r“ du — e, e —> 0+ 


and to collapse or project the wave functional right before T M , namely, collapse on a time slice almost overlapping 
the future light cone of the measurement event by Alice. It is guaranteed that there exists some coordinate system 
having such a spacelike hypersurface that intersects the worldline of Alice at r A (Tj) and the worldline of Bob at 
t b = T° dv — e in a relativistic system. 

If we further assume that the mutual influences are nonsingular and Bob performs the local operation right after 
the classical information from Alice is received, namely, at t b = Tf <to + e with e — i 0+, then the continuous evolution 
of the reduced state of detector B from t b (Tm ) = T° dv — e to T% dv + e is negligible. In this case we can calculate the 
best averaged FiQT using Eq. (17). 
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[40] The Unruh effect is both an environmental as well as a kinematic effect, the former referring to the interaction of a detector 
with a quantum field and the latter referring to its uniform acceleration (UA). Both the quantum held and the acceleration 
aspects can be referred to as relativistic, as the Unruh effect is often referred to. However, in this paper, for the sake of 
conceptual clarity we will reserve the word “relativistic” to refer to special relativistic effects between inertial frames, as 
depicted in the previous subsection plus noninertial effects as in the twin trajectories. The Unruh effect will be referred to 
specifically for the case of UA in the Minkowski vacuum, noninertial as it certainly is, in which thermality in the detector 
persists in the duration of UA. 

[41] We say a held state “collapses on a hypersurface” if the quantum state of the held degrees of freedom defined on that 
hypersurface is collapsed by a projective measurement. 

[42] In Refs. fiHl [1291 [301 this is simply called “the hdelity of teleportation,” with the averaging understood. 

[43] The left-hand side of Eq.(58) in Ref. [18] and the corresponding expression in the statements above it should be corrected 
to 16ReJ r o_Re7 r 2 + — 7i 2 /4. 


